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We derive analytical forms for non-flow contributions from cluster correlation to two-particle 
elliptic flow (i)2{2}) measure. We estimate non-flow contribution from p — > tt + tt~ decays and find 
it is negative but not a major contributor to non-flow effect in i>2{2}. We also estimate non-flow 
contribution from the recent STAR measurement of two-particle angular correlations. 
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Azimuthal distribution of charged hadrons in non-central heavy- ion collisions is anisotropic [lj. Expressed in 
Fourier series, the leading anisotropic term is the second harmonic, called elliptic flow (i>2)- Anisotropics may stem 
from early stage hydrodynamic expansion (hydro-flow) of the bulk medium created in relativistic heavy-ion collisions, 
thus making their measurement interesting |2j. However, other mechanisms, of non-flow nature, such as jet-correlation 
^ ' and resonance decays, can play a significant role [3|. A significant effort of anisotropic flow studies at RHIC is to 
1) , investigate the magnitude of non-flow effects. 

There are several methods to measure vi [J] ; all of them are affected by non-flow with differing degrees of sensitivity. 
One method, called the event plane method, is to construct the event plane ^ep from all charged particles except 
those of interest (exploiting the very fact that particles are anisotropically distributed) and calculate «2{RP} = 
(cos 2(0 — ^ep))/(cos 2A , 0ep) for particles of interest, where (cos 2A-0ep) is the event plane resolution. This method 
is affected by non-flow in the interested particles as well as in all other particles used in the event plane construction. 
Another method, called two-particle method, is to calculate ^{^l = (cos2A</>) using all particle pairs of interest 
where A(f> is pair opening azimuthal angle. This method is affected by non-flow in the interested particles. It can be 
shown that the ^{EP} and «2{2} are approximately equal 0, [l6j]. The third method, called four-particle method, 
is to obtain t>2{4} from four-particle cumulant [5|. This method is less affected by non-flow from particle clustering 
because the number of clusters with more than four particles is reduced and because non-flow effect is diluted by 
particle multiplicity to a high power. This method is not affected by resonance decays. On the other hand, flow 
fluctuations, mainly due to initial geometry eccentricity fluctuations, have different effects on i>2{2} («2{EP}) and 
«2{4} [fUzl- The fluctuation effect is positive in «2{2} but negative in «2{4}. The difference between t>2{2} and ^{4} 
is therefore a net effect of non-flow and flow fluctuations. 

In this paper we investigate non-flow effect originating from clusters. We focus on only W2{2} for two primary 
reasons: one is that i>2{2} is related to two-particle azimuthal correlation simply by w|{2} = (cos 2A0), and the other 
is that non-flow effect in i>2{2} comes from only particles used in the correlation so it is relatively straightforward to 
disentangle. 



II. NON-FLOW EFFECT FROM CLUSTER CORRELATIONS 



Suppose an event is composed of particles from hydro-medium and clusters of various types (such as minijets and 
resonance decays). Particle pairs can be decomposed into four sources: 

• particle pairs from hydro- medium (B), 

• particle pairs from same cluster (C), 

• particle pairs between hydro-medium and clusters {X), and 

• particle pairs between clusters (Y). 

The total sum of the cosines of pair opening angles is 



^cos2A<^ =B+ ^2 C+ ^2 2X+ y > 

i^j kGcluster k^cluster (k\^k,2)€zcluster 



(1) 
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^ cos 2Afaj , 




(i^j)^hydro 


C = 


^ cos 2 Afaj , 
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X = 


cos2A^, 




jGhydro 


Y = 


^ ^ COs2A0jj. 







where; 

(2) 
(3) 
(4) 
(5) 

Here i, j are particle indices, Afaj = fa — fa, and k stands for a cluster. Below we derive analytical form for each 
source. 

A. Background flow correlation 

Hydro-background particle correlation is only from hydrodynamic anisotropic flow: 

B= Y cos2A0y ; = Thy (cos 2 Afa iy ) = Vh v v\{2} hy (6) 

(i^j)£hydro 

where Thy = (Nh y {Nh y — 1)) is the number of background pairs. 

B. Particle correlation within cluster 

Particle correlation within cluster is given by 

/•27T f> '2tT p 2"7T 

C= V cos2Afaj= V d (fa)pci(fa)dfa f d (Afa,fa)dAfa f d (Afa,fa)dAfa cos2(Afa - Afa). (7) 
(i^j)ek Ja Jo Jo 

Here Afa_j — faj — fa; is the azimuthal angles of particles in cluster k relative to cluster axis fa (which can be 
defined just for convenience); fd(Afafa) is the correlation function of (daughter) particles inside cluster k relative 
to the cluster axis fa,, generally dependent of the cluster axis fa = fa — V'rp relative to the reaction plane, and 
L' fd(Afa fa)A(j) = 1; Vd(4>k) is number of (daughter) particle pairs in cluster k, generally dependent of the cluster 
axis; p c i(4>k) is the density function of cluster k relative to the reaction plane, which we will assume is given by elliptic 
flow of clusters, and p c i{fa)dfa = 1. Note, there can be many types of clusters (e.g. jet-correlation, resonance 
decays); the subscript 'd' stands for one type of clusters and we have omitted summation over all types of clusters 
from the formulism; in this work we will discuss only one type of clusters at a time. 
In general, 

C = V d (cos2Afa j ) c i (8) 

where V d is average number of pairs per cluster and (cos 2Afaj) c i is the average cosine of twice pair opening angle in 
the cluster. 

If particles inside cluster are independent of each other except all of them are correlated with the cluster axis, then 
we can factorize the correlation terms and obtain 

2 , .o„ s 2 



C = J o MMPci(fa)dfa (jf / d (A<^ fe )cos2A<^A<Aj f d (Afafa)sm2AfaiA<j>^ 

= f V d (fa)[(co S 2Afal k + ( S m2Afal k )p cl (fa)dfa. (9) 

Here (cos2A</>) ( ^ and (sin2A0)^ fc (A<f> = <f> — fa) are averages within cluster k, and are generally dependent of the 
cluster axis fa. 
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If the cluster correlation function f d (A(f),(j) k ) is symmetric about A<fi — 0, then (sin2A</>)^ = 0. Further, in 

the special case where particle correlation in clusters does not vary with cluster location <fi k , i.e., V d — const, and 
f d {A(j), 4> k ) = const., then 

C = V d (cos2Aq>) 2 cl . (10) 

C. Background-cluster correlation 

Correlation between cluster particles and hydro-medium particles is given by 

^ p2n 

X = ^ ^ COs2A0.y = / N hy (<t>k=l,2,...)phy(<t>hy)d<t>hy X 
idk j^hydro 

r 2ir p2n 

Pd((t>k)d4>k N d {4>k)fd{A4>i,<i>k)dA<l)i [cos2(0i - <j> hy )} . (11) 

JO 

Here 4>hy = 4>hy — ?Arp, and phy(4>hy) is the density function of hydro particles relative to the reaction plane (i.e., 
anisotropic hydro flow); N d (4>k) is number of (daughter) particles in cluster and is generally dependent of the cluster 
axis <pk- For generality we have taken the number of hydro-medium particles Nh y {4>k=i.2...) to depend on positions of 
all clusters. Such dependence can arise in real data analysis, such as jet-correlation analysis, from interplay between 
centrality cut and biases due to selection of specific clusters. Rewriting 4>i — 4>hy — A(f>i + 4>k — 4>hy, we have 

p2ir 

X = I N hy (0 k = la ,...)phy(4>hy)cOS2(f> h yd(t) h y X 
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2tt />2-k 

Pd(4>k)d4> k / N d ($ k )f d (A<f> i ,4> k )cos2(A4> i + 4> k )dA<j> i . (12) 



10 Jo 

f2-n r2ir 



Here we have used / dcf> k p c i(cf> k ) / dA<fr k f d (A<fr k , (f> k ) sin2(A^fc + <f> k ) = because of symmetries f d (A(j) k ,^ k ) 
J° Jo 

fd{—Atfi k , — 4> k ) and p c i{4>k) — Pci(—4>k)- Note, due to elliptic flow of clusters, cluster particles acquire elliptic flow 



^ i>2tt r2TT 



V2.d = (cos2(0 - Vrp)) = tt / pd(<}>k)d<t>k N d ((j)k)f d {A^i,(j)k)cos2(A(j) i +(j) k )dA(j)i. (13) 

Jo JO 



Using the notation in Eq. (fT3|) . we have 

X = N h yN d v 2 ,hyV2,d = N hy N d v 2 {2} hy v 2 {2}d. (14) 

Here the product of the V2's includes flow fluctuation, and equals to the product of two-particle V2 J s. This is because 
«2{2} of hydro-particles and cluster particles contain only fluctuation; non-flow does not exist between hydro-particles, 
nor between particles from different clusters. (Note, two 'clusters' can originate from a common ancestor, such as jet 
fragmentation into two p mesons which in turn decay into two pairs of pions. In our formulism, such 'clusters' are 
considered to be parts of a single cluster rather than two /9-decay clusters.) 

Again, in the special case where particle correlation in clusters does not vary with cluster location (f> k (Nd = const., 
N hy = const., and f d (A<p,<f> k ) = f d (A</))), Eq. ((T3j) becomes 

p2w p2n 

v 2 ,d = I p c i(4>k)d<t> k cos2cf> k I f d (A(j) i )cos2A<j) i dA(l) i = v2,ci(cos2A<j)) c i, (15) 

and we have 



o Jo 



X = N hy N d v 2 {2} hy v 2 {2} c i{coa2A4>) cl , 



(16) 
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D. Particle correlation between clusters 

Correlation between particles from different clusters is given by 



y = ^ ^ cos2A(fe = / p c i(^fei)#fe! / iVd(^fc 1 )/rf(A0 4 ,0 fcl )dA(/) 4 x 

Pci{4>k 2 )d4>k 2 N d (4> k2 )f d (A^ j ,4> k2 )dA<j) j [cos2(^ i - c/jj)}, (17) 



where &i and k 2 stand for two clusters. Rewriting 4>i — <f>j = A(f>i + (j> kl — A<j)j — 4>k 2 , we obtain 

Y = N 2 d vl A = N 2 v 2 2 {2} d , (18) 

where v 2 , d is given by Eq. (Til?]) . Again the cluster particle elliptic flow squared in Eq. (|18p contains flow fluctuation. 
In the special case where particle correlation in clusters does not vary with cluster location we have 

Y = A>f{2} ci <cos2A^. (19) 

E. Summary of non-flow effect from cluster correlations 

To summarize, let us now obtain the relationship between two-particle elliptic flow v 2 {2} that is affected by non-flow, 
and the real hydro-type two-particle elliptic flow V2{2}f ly . Assuming Poisson statistics, Eq. ([1} gives 

N 2 v 2 2 {2} = N 2 y vi{2} hy + N cl N 2 (cos2A<fe) d + 2N hy N cl N d v 2 {2} hy v 2 {2} d + N cl (N cl - l)N 2 v 2 2 {2} d 

= {N hy v 2 {2} hy + N cl N d v 2 {2} d ) 2 + N cl N 2 «cos2A<fe) d - v 2 {2} d ) , (20) 

where N = N) ly + N c iN d , N c i is average number of clusters, and we have taken distributions of total multiplicity and 
number of particles per cluster to be Poisson, so that (N(N — 1)) = N 2 and (N d (N d — 1)) = Aj . We have taken the 
number of cluster pairs to be N c i(N c i — 1) (i.e., not Poisson) so that the total number of pairs adds up to A 2 
Rearranging, we have 

= (^v 2 {2} hy + ^ V2 {2} d ^ + «cos2A^) cJ - v 2 {2} d ) . (21) 

For many cluster types, Eq. (|2"Tj) is generalized to 

«2{2}= (^v 2 {2} hy + ^^v 2 {2} d ) +^^0-((cos2A^) cl -v 2 {2} d ). (22) 

\ cl / cl 

We shall focus on the special case where all clusters are of the same type and particle correlation in clusters does 
not vary with cluster axis relative to the reaction plane. Using Eq. (fT5)) , Eq. (l2~Tj) becomes 

«2 2 {2} = (j^v 2 {2} hy + ^^« 3 {2} cI <cos2A^) cI ) + ^0- «cos2A<^)d - « 2 2 {2} c *<cos2A<^) . (23) 

Eq. (|2"3")l can be rewritten into 

v 2 {2} = v l{2} hy + 2^^^v 2 {2} hy (v 2 {2} cl (cos2A^ cl -v 2 {2} hy ) + 

^^) 2 (« 2 2 {2} ci (cos2A0) 2 i - v 2 {2} hy ) + ^0- «coB2A^)d " «2{2}ci<cos2A^) . (24) 

The second term on the r.h.s. is non-flow (beyond hydro-flow) due to correlation between hydro-particles and clus- 
ter particles in excess of that between two hydro-particles, and the third term is that due to correlation between 
particles from different clusters. These non-flow contributions, which are beyond hydro-flow, can be positive or 
negative, depending on the relative magnitudes of background particle flow and cluster flow diluted by particle 
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spread inside cluster. The non-flow contributions are positive when t>2{2} c ;(cos 2A0) C ; > V2{2}h y and negative when 
W2{2} c ;(cos2A0) c ; < V2{2}h y . This can be easily understood because if V2{2} c i = V2{2}h y , then the angular smearing 
of particles inside each cluster, (cos2A0) c ;, makes the angular variation of cluster particles less than that of hydro- 
particles, resulting in a negative non-flow contribution. If the net effect of cluster anisotropy and particle distribution 
inside clusters, D2{2} c ;(cos2A0) c /, equals to hydro anisotropy, then cluster particles and hydro-particles have the same 
angular variation relative to the reaction plane, resulting in zero non-flow from cross-pairs between hydro-particles 
and cluster particles and between particles from different clusters. 

The second part of the last term of Eq. (fM|) r.h.s., N( ^ d f|{2} c ;(cos 2A0) 2 ; , arises from assumptions of Poisson 
statistics, and can be safely neglected because generally w|{2} c ; << I. The first part of the last term of Eq. r.h.s., 
N % J 2 d ( cos 2 A0jj ) ci , is non-flow due to correlation between particles in the same cluster. This non-flow contribution 
can also be positive or negative. If particle emissions within clusters are independent, (cos 2A0 i:) ) c ; = (cos2A</>) 2 ; , 
then Eq. (pMf becomes 

vi{2} = v 2 2 {2} hy + 2 N ^ Ncl ^ d v 2 {2} hy ( W2 {2} c; (cos2A0) d - v 2 {2} hy ) + 

(^) 2 H{2Mcos2A«^ - vl{2} hy ) + ^0- (I - vl{2} cl ) (cos2A«^. (25) 

In this case the non-flow contribution due to particle correlation within clusters can only be positive. 

We note that the non-flow contributions from the second and third term of Eq. ([24| r.h.s. have the identical 
azimuthal shape relative to the reaction plane as that of hydro-flow, because they arise from the common correlation 
of clusters and hydro-particles to the reaction plane. As a result these non-flow contributions will unlikely be separated 
from medium hydro-flow in inclusive measurement of azimuthal correlation. To separate these two contributions, one 
needs to identify clusters and measure two-cluster azimuthal correlation. In fact, elliptic flow is often defined as 
the second harmonic of particle distribution relative to the reaction plane, t>2{RP} = (cos 2(0 — ?/>rp)) . For events 
composed of hydro-particles and clusters, we have 

z, 2 {RP} = ^(cos2(0-VRp)) % + ^^<cos2(0-^ R p)} d = ^«a{RP} hv +^^w 2 {RP} cI <coB2A0) cI . (26) 

This is analogous to the terms in the first pair of parentheses on Eq. (|23p r.h.s. except the latter contains flow 
fluctuation. The elliptic flow definition by Eq. (|26[) contains cluster contribution through angular spread of particles 
in clusters, (cos2A</>) c /, and anisotropy of the clusters themselves, v 2 {RP}ci- This raises question to comparisons 
often made between elliptic flow measurements and hydro calculations which may include flow fluctuation but does 
not include cluster correlations. 

III. ESTIMATE OF NON-FLOW FROM TWO-BODY RESONANCE DECAYS 

In relativistic heavy-ion collisions, a large fraction of final state pions come from p decays. Charged pion pairs 
from p — > 7r + 7r~ decays have intrinsic angular correlation due to decay kinematics. In this section, we estimate the 
non-flow effect from this intrinsic angular correlation on pion elliptic flow. 

Suppose a parent particle of mass M decays into two daughter particles of equal mass m with decay angle 9 (the 
angle a daughter particle makes in the parent cm. frame with the parent direction of motion). The opening angle 
between the two daughters is straightforward to calculate: 

87 2 /S 2 (l -4a 2 ) sin 2 6» 

cos2A0i2 = 1 1 —— - 5 , (27 

[ 7 2 + 4a 2 - 7 2 /3 2 (I - 4a 2 ) cos 2 8} 2 - 16 7 2 a 2 



where (3 is the parent speed, 7 = I/yl — /3 2 its Lorentz factor, and a = m/M. To obtain the opening angle projected 
onto the transverse plane, however, is a bit tedious. So we resort to MC sampling of p — > decays. We assume 

uniform rapidity distribution for p within \y\ < 2. We calculate the average cosine of twice the opening azimuthal 
angle, (cos2A0i2) p , between the two daughter pions in the lab frame if both are within rapidity |y| < I. We assume 
isotropic decay in p cm. frame 

Figure I (left panel) shows (cos 2A0i 2 ) p as a function of p transverse momentum (pr)- (cos 2A0i 2 ) p = I for px = 0, 
because the decay pions are back-to-back in the lab frame, (cos 2A(/>i 2 ) p drops quickly with pt due to closing of pion 
pairs from boosting. (cos2A0i2) p becomes negative at intermediate pt when the opening angle is narrowed towards 
7r/2. When increases further, strong boost focusing narrows the pair opening angle to asymptotically approaching 
zero. 
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FIG. 1: Left panel: (cos 2A(j>\2) P vs parent p-meson pt- Flat pr distribution is used to have uniform statistics. Right panel: 
(cos2A0i2) P vs daughter pion pr- Each decay is entered twice in this plot, at pr of each pion. The transverse distribution is 
assumed to be m-r-exponential with inverse slope T. Rapidity range \y\ < 2 is simulated for p and daughter pions are restricted 
within \y\ < 1 for both panels. Uniform decay angle distribution is used. One million decays are simulated for each set of data 
points. 



TABLE I: p — > tx + -k decay (cos2A<^i2) p versus T (p mr exponential inverse slope). 



T (GeV) 0.1 


0.2 


0.3 


0.4 


0.5 


0.7 1.0 


(cos2A0i 2 )p 0.13 


-0.13 


-0.17 


-0.14 


-0.10 


-0.03 0.05 


(cos2A<^) p 0.07 


0.18 


0.28 


0.36 


0.41 


0.49 0.55 



To calculate the average (cos2A0i 2 )„over all pr, we take exponential in mr — \/ M 2 + p T for transverse distribu- 
tion as measured in Au+Au collisions [lOj : 

( mT \ (28) 

Figure 1 (right panel) shows (cos 2A</>i2) p as a function of daughter pion pr- Each decay is entered twice in this 
plot, at pt of each pion. Table I lists the overall (cos2A0i2) p versus T. (cos2Acj)i 2 ) decreases and then increases 
with T. In the range T — 0.2 — 0.5 GeV relevant for heavy-ion collision [T(J], (cos2A</!>i2) p is negative and relatively 
constant. Emissions of the two daughter particles are not independent, so (cos 2A(/>i 2 ) p ^ (cos2A0) 2 , where Acj>i 2 
is azimuthal opening angle of the daughter pion pair and Aip is azimuthal angle of daughter pion relative parent 
direction of motion. Table I lists (cos2A0) p together with (cos2Acf>i 2 )p at various T values. 

Let x be the fraction of final state pions from p decays; the number of p's (clusters) is xN/2. We will assume 
Poisson statistics except for the number of daughters (Nk = 2) and daughter pairs (Nk(Nh — 1) = 2) per decay that 
are both fixed. Because of the non-Poisson statistics of the p-decay daughters, Eq. (j24]l cannot be readily applied. 
However, this can be easily revamped by noting that the individual non-flow contributions in Eq. ([T]) are now: 

B = {l~x) 2 N 2 v 2 {2} hyi (29) 
C = 2(cos2A</> 12 ) p , (30) 
X = 2{l-x)Nv 2 {2} hy v 2 {2} p {cos2A(f>) p , (31) 
Y = 4?; 2 {2} p (cos2A</>) 2 . (32) 

Thus Eq. |T]) becomes 

N 2 v%{2} = (l-x) 2 Af 2 vf{2} /ly -kciV(cos2A ( ^ 

(33) 

where we have taken the number of cluster pairs to be f 2 ^) • ( xN 2 1 ) so that the total number of particle pairs adds 
up to A^ 2 . Finally we obtain 



v 2 {2} = v 2 {2} hv + 2xv 2 {2} hy (v 2 {2} p (cos2A(j)} p -v 2 {2} hv ) + 
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x 2 (v 2 {2} p (cos2A0) p - v 2 {2} hy ) 2 + — ((cos2A0 12 ) p - u|{2} p (cos 2A<^) 

= [(1 - x)v 2 {2} hy + OT2 {2} p (cos2A0) p ] 2 + ^ ((cos2A0i 2 ) p - w|{2} p (cos2A0)2) . (34) 

The second and third terms of Eq. (|3~4")l r.h.s. are non-flow contributions from correlations between decay particles 
and hydro-particles and between daughter particles from different decays. Obviously, the magnitude of these non-flow 
depends on ii 2 {2} p : if w 2 {2} p (cos 2A<j)) p > V2{2}f ly , then their non-flow is positive, and otherwise, negative. The 
last term of Eq. ()34|) r.h.s. is the non-flow contribution from correlation between the two decay daughters, and is 
negative for p — > tt + tt~ decays with our chosen kinematics. Since p is short-lived, p and hydro-medium pions are in 
detailed balance at the early stage of the collision where elliptic flow is generated. It is therefore reasonable to assume 
that w 2 {2} p (cos 2A(/))p = v 2 {2}h y , i.e., p decay pions and hydro-medium pions possess the same elliptic flow and are 
indistinguishable. In this case the non-flow contribution is entirely from correlations between daughter particles of 
the same decay: 

v 2 {2] = v 2 2 {2} hy + |p ((cos2A0 12 ) p - « 2 2 {2} p (cos2A0)2) » v\ hy + ^(cos2A0 12 ) p . (35) 

Here we have made assumption v 2 {2} p << 1. It is worth to note that only the final-stage p decays are relevant 
because pions from decays at early times rescatter and lose their intrinsic angular correlation. Those decay pions will 
eventually become part of the hydrodynamic medium after rescattering. The fraction of p-decay pions x, therefore, 
should refer to only those final-stage p decays that are experimentally measured. 

Table II lists our estimates of non-flow contributions from p — > tt + tt~ decays for various centrality bins together 
with V2 measurements [111 ]. The p ttit exponential inverse slope is taken to be T — 0.3 GeV as measured in peripheral 
Au+Au collisions [l(j and consequently (cos2A0i 2 ) p = —0.17; other T values reasonable for heavy-ion collision do 
not alter (cos 2A<fti2)n significantly. The fraction of decay pions is taken to be x = 40% as measured in peripheral 
Au+Au collisions [10| : its centrality dependence is weak as the chemical freeze-out temperature is measured to be 
independent of centrality [l2| and losses due to final state rescattering are likely offset by regeneration processes [l3[ . 
The listed V2{2}h y values are calculated from the measured f 2 {2} assuming the only non-flow contribution is from 
p — > 7r + 7r~ decays and w 2 {2} p (cos 2A<p) p — V2{2} hy Magnitude of non-flow contribution is calculated by 1 — ^rfe 2 , 
and is typically a few percent for most centrality bins. This is small (and with wrong sign) compared to the difference 
between the w 2 {2} and u 2 {4} measurements (see 1 — ^| 2 j in Table II). Note that this non-flow contribution is only 
from correlation between daughter particles from the same decay, whilst contributions from correlations between 
particles from different decays and between decay particles and medium hydro-particles are arguably small (i.e., 
« 2 {2} p (cos2A0) p = V2{2}h y ) and are neglected here. Even not small, those contributions should be already included 
in the w 2 {4} measurement. Therefore, the large difference between the measured v 2 {2} and u 2 {4} suggests that 
resonance decays are not a major contributor to non-flow at RHIC. The difference must then come from other major 
non-flow effects and/or flow fluctuations resulting from initial eccentricity fluctuations. 



IV. NON-FLOW FROM "MINUET" CORRELATIONS 

Strong angular correlations in (A77, A0) between soft particles have been observed by STAR [l5j]. They reveal 
in part characteristics of jet-like correlations from small energy partons (the so-called minijets), such as the narrow 
A(f> correlations. On the other hand, the correlations are measured to be extraordinarily wide in Arj. The long Arj 
correlation ridge may be due entirely to the medium, such as fluctuations of color flux tubes in the initial stage of the 
collisions. The ridge particles are then focused into the narrow A<j) window by the radial flow boost, a result of the 
medium's hydrodynamic expansion. This would suggest that the measured two-particle angular correlations may not 
necessarily be due to jet-like correlations, but hydrodynamic expansion of the medium. Even if the measured angular 
correlations are indeed due to jet-like correlations, those correlations must have been influenced significantly by the 
medium hydrodynamic flow. It is, therefore, difficult to separate flow and non-flow. 

The measured angular correlation raw data (i.e., pair density per charged hadron) were decomposed into two 
components: elliptic flow and minijet correlations. The latter is assumed to consist of a two-dimensional Gaussian in 
(A77, A(f>) and a dipole moment in A<p: 



f(Acj), A77) oc A A ^ cos A0 + Ai cxp ( ~ ) ■ (36) 



At,/ 

We note, however, that this decomposition is strongly model dependent. If the assumed functional form of Eq. (|3"o]l 
is incorrect for minijet correlation, then the extracted minijet correlation will contain flow contributions, as discussed 
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TABLE II: Estimates of non-flow contributions from p — > tt + tv~ decays in Au+Au collisions at 200 GeV as a function of centrality 
together with two-particle «2{2} and four-particle i>2{4} measurements [jj. The dN^/dri results are from Ref. The p 

transverse distribution is taken to be ttlt exponential with an inverse slope of 0.3 GeV [101 ] (and consequently (cos2A0i2) p = 
—0.17), and the fraction of pions from p —¥ tt + tt~ decays is taken to be 40% [lb] ], both independent of centrality. It is also 
assumed that decay pions and bulk hydro-medium pions possess the same elliptic flow, i.e., V2{2} p (cos 2A(f>) p = V2{2}h y , and are 
indistinguishable. V2{2}h y is calculated from the measured V2{2} assuming the only non-flow contribution is from p — > -k + -k~ 

decays. Relative non-flow contribution from resonance decays is calculated as 1 — -fjjg-p-, to be compared to non-flow effect 

indicated by 1 — V2 \t] from measurement. 

j) 2 {2} 
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0-5% 
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2.41% 




-0.70% 


2.51% 


-4.15% 




5-10% 


558 


3.55% 


2.53% 


-0.78% 


3.63% 


-2.39% 


28.73% 


10-20% 


421 


4.97% 


4.27% 


-0.90% 


5.05% 


-1.62% 


14.08% 


20-30% 


287 


6.42% 


5.66% 


-1.09% 


6.51% 


-1.43% 


11.84% 


30-40% 


195 


7.29% 


6.33% 


-1.32% 


7.41% 


-1.63% 


13.17% 


40-50% 


126 


7.64% 


6.43% 


-1.64% 


7.81% 


-2.29% 


15.84% 


50-60% 


78 


7.59% 


6.18% 


-2.09% 


7.87% 


-3.71% 


18.58% 


60-70% 


45 


7.25% 


5.68% 


-2.75% 


7.75% 


-6.95% 


21.66% 


70-80% 


22 


6.88% 




-3.93% 


7.92% 


-15.17% 





above. Therefore, the name "minijet" here should be taken merely as a label to refer to all the correlation structures 
except elliptic flow under the particular assumption of Eq. (|36[) . Keeping in mind the model dependence, we now 
proceed to estimate the magnitude of non-flow contributions from the decomposed minijet correlations. 

We shall assume minijet correlation does not vary with minijet direction relative to the reaction plane and particle 
emission within minijets is independent, and use Eq. (|25p to estimate the magnitude of minijet contribution to non- 
flow from STAR's measurement [lf|. From Eq. the average pair opening angle can be obtained by (cos2A0) = 
J Q * f(A(f>) cos 2A0g?A</> and equals to (cos2A0)^ due to independent particle emission. Number of minijet particle- 
pairs per hadron (V) can be obtained by integrating f(A(f>, Ary) within the STAR acceptance (2 units of pseudo- 
rapidity 77). Total number of minijet particle pairs within acceptance is then N c iN% = VN where TV = 2dN c ±/ drj. 
To obtain the number of minijet particles per hadron N^N^/N , we note that the measured number of particle pairs 
per minijet is VN/N c i where N c i is number of clusters contributing to the measured signal, and in turn number of 
fragments per hadron is N cl N d /N = N c iy/VN/N c i/N = ^VN c i/N. Number of minijets per unit of pseudo-rapidity 
should be just that in pp scaled by the number of binary collisions, 0.013A^bm- However, minijets outside acceptance 
also contribute to the measured signal, thus N c i > 0.013A^bi n ; perhaps N c i — 0.013iVbin X k where n ~ 2 given the 
measured broad 77-width of minijet correlation. (Note, minijets inside acceptance always contribute to the measured 
signal, although some of the minijet fragments will leak out of acceptance, which makes the measured fragment 
multiplicity per minijet smaller.) Number of fragments per hadron is therefore N c iNd/N > ^/0.013ATbinV7 {dN c h j drj) , 
good to a factor of i/k. The cluster size (fragment multiplicity per minijet) measured inside acceptance is in turn 
N d — V/(N c iN d /N) < ^yV(dN c h/drj)/(0.013Nb in ), again good to a factor of ^/n. In the estimates below, we assume 

K = 1. 

In order to estimate non-flow effects arising from cross-talk between cluster correlation and cluster flow, terms X 
and Y in Eq. (|4]) and ([5]), respectively, we need to know the elliptic flow of clusters. This is not measured. We may 
estimate its magnitude at least in two ways, and they can give very different results, (i) Minijet cluster size (not for 
clusters measured in acceptance, but those in full space) is ^JV{dN c h/drj) / (0.013iVbin) which is roughly 10 charged 
hadrons in central Au+Au collisions, suggesting the parent px ~ 6 GeV/c. Those initial partons may have large v-i 
(taking on the saturated value at high px), which is x3-4 higher than the average, (ii) If clusters result from initial 
state (semi-)hard parton scatterings, then their distribution should be isotropic, vi — 0. Subsequent jet quenching 
results in cluster size varying with the reaction plane, larger out-of-plane than in-plane. This yields effectively a 
finite cluster v-i (weighted by the number of particle pairs per cluster) which can be negative. Thus, depending on 
physics scenario, cluster V2 can be positive or negative. However, non-flow effects from cross-talk terms are generally 
small, significantly smaller than that from particle correlations within the same cluster, i.e., C in Eq. ([3]). Thus, 
we shall neglect non-flow effects from cross-talk in the following estimate. Again, it is worth to note that non-flow 
due to minijet-hydro and minijet-minijet correlations have identical angular shape as elliptic flow, so they cannot be 
distinguished in experiment and thus are included in essentially all elliptic flow measurements. 

We estimate non-flow effect for top 5% central Au+Au collisions at 200 GeV. The measured minijet amplitude is 
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Ai « 0.63, T^-width is sa 2.1, and total extrapolated volume is 27r^4iCTA7)CrA0 ~ 5.3 [1 51 ] - From these measurements 
we calculate the minijet </>-width to be <7a<£ ~ 0.64, the minijet azimuthal angle spread to be (cos 2Acf>) — (cos 2A0)^ « 
0.44, and the minijet signal volume within STAR acceptance to be V ~ 1.9. The number of minijet fragments is 
N cl N d /N w ^0.013A^binV/(diV ch /dr/) w 0.19 where dN ch /dr) = 691 and iV bin = 1012 [3]. Non-flow contribution from 
minijet correlation is given by ^ C = (F/iV)(cos 2A0) « 0.025. Comparison to the measured «2{2} = 0.024 by the 
two-particle method 5] suggests that non-flow from correlation between particles from the same minijet cluster is the 
dominate contributor to v 2 {2}. 

We note, however, that the above estimate and conclusion are based on the assumption that Eq. (|3l)f is the proper- 
functional form for mini jet correlation so that the decomposition of the measured two-particle correlation into minijet 
and elliptic flow in Ref. (T5J is correct. It is possible that our result may be an overestimate of non-flow contributions 
because, as discussed earlier, the minijet correlation may contain flow contributions and is most likely affected by 
hydrodynamic how. In fact, ample evidence indicates that flow fluctuation effects are significant in central Au+Au 
collisions 6, 7]. The fact that our estimated non-flow contribution from minijet correlations dominates the measured 
t>2{2} may suggest that the decomposition of elliptic flow and minijet correlations by Eq. (|36[) is improper resulting 
in an overestimate of the non-flow. 



V. SUMMARY 

We have derived analytical forms for non-flow contributions from cluster correlations to two-particle elliptic flow 
measure «2{2}. We estimate non-flow contribution from two-body p — > tt + h~ decays. With transverse distribution in 
accordance with measurement, non-flow contribution from p decays is negative, contradictory to common perception. 
The magnitude of the non-flow contribution is small, on the order of a few percent for most centrality bins of 
Au+Au collisions. The large difference between two- and four-particle elliptic flow measurements cannot be due to 
resonance decays. The likely sources for the difference are eccentricity fluctuations and other non-flow effects such as 
jet correlations. 

STAR has decomposed the two-particle (A?y, A<j>) correlation into a quadrupole component t>2{2D} and the so- 
called minijet correlations assuming the functional form of Eq. (|3"6| for the latter [r| . We outlined a procedure to 
estimate the non-flow contributions from the decomposed minijet correlations to elliptic flow «2{2} measurement by 
the two-particle method. We estimated the magnitude of the non-flow contributions in central Au+Au collisions, and 
found the non-flow contributions from the decomposed minijet correlations to be predominant. This may suggest, 
because flow fluctuation effects are expected to be sizeable in central Au+Au collisions, that the assumed functional 
form for minijet correlations used in the decomposition may be improper. 

Nevertheless, given the decomposed minijet cluster correlation, we found the main non-flow contribution is from 
correlations between particles within the same cluster. The magnitudes of non-flow contributions due to correlations 
between cluster particle and hydro-particle and between particles from different clusters are generally small, depending 
on the relative vi of hydro-particles and clusters themselves. Those non-flow effects from cross-talk of particles are 
included in all available elliptic flow measurements. Thus comparisons between measurements and hydro calculations 
should be taken with caution. 
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